Multi-Phase Patterns in Periodically Forced Oscillatory Systems 
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Periodic forcing of an oscillatory system produces frequency locking bands within which the system 
frequency is rationally related to the forcing frequency. We study extended oscillatory systems that 
respond to uniform periodic forcing at one quarter of the forcing frequency (the 4:1 resonance) . These 
systems possess four coexisting stable states, corresponding to uniform oscillations with successive 
phase shifts of 7r/2. Using an amplitude equation approach near a Hopf bifurcation to uniform 
oscillations, we study front solutions connecting different phase states. These solutions divide into 
two groups: 7r-fronts separating states with a phase shift of vr and 7r/2- fronts separating states with a 
phase shift of 7r/2. We find a new type of front instability where a stationary 7r-front "decomposes" 
into a pair of traveling 7r/2-fronts as the forcing strength is decreased. The instability is degenerate 
for an amplitude equation with cubic nonlinearities. At the instability point a continuous family of 
pair solutions exists, consisting of 7r/2-fronts separated by distances ranging from zero to infinity. 
Quintic nonlinearities lift the degeneracy at the instability point but do not change the basic nature 
of the instability. We conjecture the existence of similar instabilities in higher 2n:l resonances 
(n = 3,4,..) where stationary tt- fronts decompose into n traveling vr/n- fronts. The instabilities 
designate transitions from stationary two-phase patterns to traveling 2n-phase patterns. As an 
example, we demonstrate with a numerical solution the collapse of a four-phase spiral wave into a 
stationary two-phase pattern as the forcing strength within the 4:1 resonance is increased. 



I. INTRODUCTION 

Periodic forcing of an oscillatory system produces a 
multiplicity of uniform stable phase states. The simplest 
situation arises within the 2:1 frequency locking band 
where the system oscillates at one half of the forcing fre- 
quency. In that case "two-phase" patterns appear, in- 
volving alternating domains that oscillate with a phase 
shift of TT ^-|^. The boundaries between nearby domains, 
hereafter tt- fronts, may undergo a parity breaking bifur- 
cation rendering a stationary front unstable and giving 
rise to a pair of counterpropagating fronts M. This in- 
stability, the so called "nonequilibrium Ising-Bloch bifur- 
cation" (or NIB), designates a transition from standing 
two-phase patterns to traveling two-phase patterns P-R ■ 
The instability is demonstrated in Fig. H as a grey-scale 
map in the space-time plane. Recent experiments on 
a photo-sensitive Belousov-Zhabotinsky (BZ) reaction, 
periodically illuminated, have also revealed a transition 
to labyrinthine patterns within the 2:1 band, suggest- 
ing the possible existence of a transverse instability of 
TT-fronts 1). 

The situation becomes more complicated within the 
4:1 band which has four stable phase states shifted by 
7r/2 with respect to one another M. In addition to n- 



fronts there also exist 7r/2-fronts separating oscillating 
domains with a phase shift of 7r/2. The multiplicity of 
front solutions increases with the order of the band. The 
6:1 band has three types of fronts: 7r-fronts, 27r/3-fronts 
and 7r/3-fronts. The 8:1 band has four types of fronts (tt, 
37r/4, 7r/2 and 7r/4) and so on. In addition to adding new 
types of fronts as the band order is increased the number 
of front solutions of a given type also increases. 

In this paper we report on a new instability of tt- fronts, 
occurring within the 4:1 band. Upon decreasing the forc- 
ing strength a stationary 7r-front loses stability and de- 
composes into a pair of traveling 7r/2 fronts. The insta- 
bility is demonstrated in Fig. 0. The decomposition into 
a pair of traveling 7r/2-fronts is accompanied by the ap- 
pearance of an intermediate (grey) domain whose phase 
of oscillation is shifted by tt/2 with respect to the adja- 
cent white and black domains. Like the NIB bifurcation 
the TT- front instability within the 4:1 band designates a 
transition from stationary patterns to traveling waves. 
The significant difference is that the two-phase station- 
ary patterns give place to traveling four-phase patterns. 
This feature of the 4:1 resonance is related to a peculiar 
prop erty of the vr-front instability to be discussed in Sec- 
tion III, The TT- front decomposition instability appears 
to exist in higher 2n:l bands as well. We analyze in de- 



tail the 4:1 resonance case and bring numerical evidence 
for the existence of this type of instability in the 6:1 and 
8:1 resonances. A brief account of some of the results to 
be reported here has appeared in Ref. pO|. 
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FIG. 1. The NIB bifurcation in the 2:1 resonance: 

space-time plots showing an unstable stationary vr-front 
(Ising) evolving into left (a) and right (b) traveling vr-fronts 
(Bloch) beyond the NIB bifurcation. 
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FIG. 2. The decomposition instability in the 4:1 resonance: 
Space-time plots (solutions of Eqn. iMj showing the decompo- 
sition of an unstable 7r-front into a pair of 7r/2-fronts traveling 
to the left (a) or to the right (b). The pairs of 7r/2-fronts en- 
close grey colored domains whose oscillation phases are shifted 
by 7r/2 with respect to the black and white domains. Param- 
eters in Eqn. (|I|): ^ = 1.0, f = 0.02, 74 = 0.3. 

We consider an extended system that is close to a Hopf 
bifurcation and externally forced with a frequency about 
four times larger than the Hopf frequency. The set of 
dynamical fields u describing the spatio-temporal state of 
the system (e.g. set of concentrations in the BZ reaction) 
can be written as u = UoAexp {i-f't) + c-c + . . ., where 
Uq is constant, A is a slowly varying complex amplitude, 
ujf is the forcing frequency and the ellipses denote smaller 
contributions. The equation for the amplitude A can be 
written in the following standard form (after rescaling 
and shifting arg ^ by a constant phase) ||ll|-p^ : 



Ar^{n + iv)A + (1 + ia)A,_^ - (1 - ili)\A\'^A 
+ 74A*' , 



(1) 



where the subscripts r and z denote partial derivatives 
with respect to time and space, and all the parameters 
are real. The proximity to the Hopf bifurcation implies 
H « 1. We will also be using the following form of 



Eqn. (||) obtained by rescaling time space and amplitude 
as t = fj,T, X = yj^/2z and B = A/^JJl: 

Bt = {l + iuo)B + \{l + ia)B^. - (1 - ^P)\B\''B 

+ liB*\ (2) 



where vq = v j [i. 

II. FRONT SOLUTIONS 

We first study the gradient version of Eqn. 
is obtained by setting vq — a = (i = ^: 

Bt=B + ]^B,,,-\B\''B + -nB*\ 



which 



(3) 



Eqn. (pi) has four stable phase states for < 74 < 1 shown 
by soUd circles in Fig. ||: B±i = ±A and B±i = ±iA, 
where A = 1/-^1 — 74. Front solutions connecting pairs 
of these states divide into two groups, 7r-fronts and 7r/2- 
fronts. The 7r-fronts, shown in Fig. y as solid lines, are 
given by 



i3_i^+i = B^i tanh x . 
B^i^^i = B^i tanh a; . 



(4) 



The 7r/2-fronts are shown in Fig. || by the dashed curves. 
For the particular parameter value 74 = 1/3 they have 
the simple forms 
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FIG. 3. Phase portrait in the complex plane of solutions 
to Eqn. (pi). The dots represent the 4 spatially uniform phase 
locked solutions. The solid lines are the tt front solutions and 
the dashed lines are the 7r/2 fronts. The thin lines in the circle 
are the phase portrait of showing the collapse of a vr front into 
two 7r/2 fronts. 
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(^0 + \\^P) < Bo > +\P < BiB', > 



-Mo + P), 



(7) 



where A = •y/3/2. A perturbation analysis around 74 = 
1/3 shows that the expression (0) for the speed remains 
valid for small deviations of 74 from 1/3. 



B+i- 



= ly2 

= 1^2 



[1 + i — (1 — i) tanh x\ , 
[1 — i + (1 + z) tanh x\ , 



(5) 



Additional front solutions follow from the invariance of 
Eqn. (0) under reflection, x — > ~x. For example, the 
symmetric counterparts of B^i^^i{x) and B^i^^i{x) 
are B+i^+i{x) = B+i^+i{-~x) and B^i^+i{x) = 
B+i^^i[-x). 

Consider now the nongradient system (0). The main 
effect of the nongradient terms is to make the 7r/2-fronts 
traveling. The nongradient terms have no effect on the 
TT-fronts which remain stationary. To see this we assume 
a traveling solution B{x — ct) of Eqn. (0) and project this 
equation on the translational mode B' . For tt- fronts we 
obtain 



c < B'o >= 



Bo(z) = A tanh z , 



(6) 



implying c — (the brackets denote integration over the 
whole line). For 7r/2-fronts with 74 — 1/3 we find 



III. A n-FRONT INSTABILITY 

The TT-fronts (Q) are similar to the Ising front in the 
2:1 band but as we will see shortly the instability they 
undergo is not a pitchfork bifurcation like the NIB. It is 
rather a degenerate instability leading to asymptotic so- 
lutions that are not smooth continuations of the unstable 
stationary 7r-fronts in a sense to be made clear in the fol- 
lowing. A stability analysis of the 7r-fronts indicates that 
they lose stability at 74 = 1/3. To analyze the instability 
we study Eqn. (0) near that critical value. 

A. Gradient system 

We begin with the gradient version (g). Introducing 
the new variables 

[/ = 5R(B) + 5(B) V = n{B)~Q{B), (8) 

we rewrite Eqn. (0) as 

C/t-C/+^C/..-^t/'-^(C/'-3y2)c/, (9a) 

Vt = V + ^V,^-^V^-^{V^-3U^)V, (9b) 



where 



d = 74 - 1/3 , 



At the instabihty point, 74 — 1/3, the two equations 
decouple (since d — 0) and admit solutions of the form 



U — aiBQ{x — xi) . 

V = (T2Bq{x~ X2) ; 

- tanha;, cri^2 



(10) 



where -60(2;) = \/| tanha;, cri^2 = ±1, and xi and X2 are 
arbitrary constants. An intuitive understanding of this 
family of solutions can be obtained by expressing these 
solutions back in terms of the complex amplitude B. For 
(Ti = —(72 = 1 for example, the solution (HG) is equivalent 
to 

B{x]Xi,X2) == i?-i^+i(x- xi) + B+i^+i{x - X2) - X. 

When |x2 — xil -^ cx) this form approaches a pair of 
isolated 7r/2-fronts: 



and 



B K B^ 



B K B. 



^+i{x-xi), 



i{x ~ X2) 



Xi , 



X2 ■ 



When X2 — xi = it reduces to the 7r-front B^i^^i. 
Defining a "center of mass" coordinate, C, and an order 
parameter, x, by 



C = -^(xi +X2). 



X^^[X2 



{X2 - Xi), 



the one-parameter family of solutions, {B{x;(^,x) I X G 
R}, where B{x; C, x) = B{x; xi, X2), represents 7r/2-front 
pairs with distances, 2%, ranging from zero to infinity. 

For I74 — 1/3| = |d| <C 1, the weak coupling between 
the two equations ( |9a| ) and (9b) induces slow drift along 
the solution family B{x; xi,X2). We now write a pair 
solution as 



U = (7iBo[x — xi{t)] + u , 
V = a2Bo[x- X2it)] +v, 



(11) 



where u and v are corrections of order d. Inserting these 
forms in Eqns. m) we obtain 



Hiu — aiXiBQ{x — xi) 



(12) 



-dai[BQ{x - xi) - Wq{x - X2)]Bo{x - xi) , 



TI.2V = CF2X2B'q{x - X2) 



(13) 



- -jdcr2[Bl{x - X2) - 3i3o(a;- xi)\Bq{x - X2) , 

where Hi, 2 = —^— ^-§^ + '^Bq{x — xi^2)- Projecting the 
right hand side of Eqn. ( |l2|) onto Bq{x — xi), the zero 
eigenmode of 7i{ — Hi, and setting to zero we obtain 



27 



dx 



X tanh(a; — xi) sech^(a; — xi) tanh^(a; — 2:2) , (14) 
A similar solvability condition for (|3|) leads to 

27 



X2 



16 



dx 



X tanh(a; — 2:2) sech (x — X2) tanh {x ~ xi) . (15) 
Expressing these equations in terms of C and x we find 






(16) 
(17) 



where 

/oo 
dz tanhz sech^z tanh^(z + 2x). (18) 
-00 

Evaluation of the integral in (08) yields 

J(x) = I{a) = 6(a'i - a-^) + (1 - 3a-^)G{a) , 
G(a) = (l-a-2)ln(i±^), 

where a == tanh2x. Note that Eqns. { {Wj) and ( [l7|) are 
valid to all orders in x and to linear order around 74 = 

1/3- ^ 

The equation for the order parameter ( |l7|) can be writ- 
ten in the gradient form 



dV 
dx' 



27 
^-16^ 



J{z)dz . 



(19) 



Fig. ^ shows the potential V{x) for d > (74 > 1/3) 
and d < 0. There is only one extremum point, x = 0, 
of V. For d > it is a minimum and x converges to 
zero. Pairs of 7r/2-fronts with arbitrary initial separa- 
tion, X2 — xi, attract one another and eventually collapse 
to a single 7r-front {xi = X2 or x = 0). In practice, 
the collapse process is noticeable only for relatively small 
separations. For d < the extremum point, x = 0, is a 
maximum and x diverges to ±00. A 7r-front decomposes 
into a pair of 7r/2-fronts which repel one another. This 
process is shown in Fig. 0for the nongradient system (1^). 
In the gradient case both it and7r/2-fronts are stationary 
(in the absence of interactions). Since the potential V{x) 
becomes practically flat at finite x values, the pair of 7r/2- 
fronts do not seem to depart from one another at long 
times. Fig. shows the decomposition process of a tt- 
front in the complex B plane. Starting with the _B_i_,+i 
TT- front, represented by the thick solid phase portrait, the 
time evolution (thin solid phase portraits) is toward the 
fixed point B^i and the dashed phase portraits represent- 
ing the pair of 7r/2-fronts i?+i^+i and i?+i^_i. Because 
of the parity symmetry x ^ ~Xi an appropriate pertur- 
bation of the initial i3_i_+-f 1 7r-front could have led the 
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FIG. 4. The potential V{x)- (a) For d > the extremum 
at X = is a minimum and and x converges to 0. (b) For 
d < the extremum is a maximum and x diverges to ±cx3. 



dynamics toward the pair i?+i^_i and _B_i^_i. Notice 
that for d = 0,(' = 0,x = Oi ^^^d we recover the two- 
parameter family of pair solutions B{x; C, x) with arbi- 
trary C and X- This degeneracy of solutions at d = is 
lifted by higher order terms in the amplitude equation p) 
as will be discussed in Section [II C below. 



-dcri[i?o(x- xi) - 3Bq{x - X2)]Bo{x - xi) 



+ vo(T2Bo{x - X2) + -ao-2-Bo (a; 



X2) 



-/3(72 [Bl{x - xi) + Bl{x - X2)\Bo{x - X2) , (21) 



H2V = a2X2BQ{x 



X2) 



-da2[Bl{x - X2) - iBl{x - xi)]So(a 



X2) 



vocriBo{x - xi) - -aaiB'^ix 



Xi) 



- -Pai[Bl{x - xi) + B^ai^ - X2)]Boix - xi) . (22) 

Solvability conditions lead to equations for x and C. The 
equation for x remains unchanged. That is, Eqn. (17) 
is valid for the nongradient equation (ph as well. The 
equation for C becomes 



aifJaC = i^oF^ix) + aF^ix) + PPpix) 



(23) 



where 



F, 

Fa, 

Fp 



--G{a) 



9-2 



a-^G{a). 



Notice that F^, Fa, and Fjj are odd functions of x and do 
not vanish when d = Q. When |x| — > cxd the right hand 
side of (|23| ) converges to |(i'o -I- /?), the speed of a 7r/2- 
front solution of Eqn. (g). The odd symmetries of Fy, 
Fa, and Fp imply that the x — solution (representing 
a TT-front) remains stationary [C, = 0) in the nongradient 
case as well, and that the two pairs of 7r/2-fronts x ~ ±00 
propagate in opposite directions. 



B. Nongradient system 

The results described above can easily be extended to 
the nongradient system (g) for small a, /3 and Vq. The 
equations for U and V are 



Ut^U+hj^^ 



3 2^ ^ 






voU-^-U^..-^{U^ + V^)U. 



a, 
— ( 
2 



(20) 



Assuming d, a, (3 and vq are of the same order of mag- 
nitude we write a solution in the form (O), insert in 
Eqs. (pq) and obtain 

Tiiu — aiXiB^ix — xi) 



C. The effect of higher order terms 

According to Eqn. ( [l7| ) the asymptotic solutions just 
below 74 = 1/3, the 7r/2-front pairs as \x\ — > 00, are 
not smooth continuations of the stationary 7r-front at 
74 = 1/3 (the X — solution). This abrupt nature of 
the instability is related to a degeneracy of solutions at 
74 = 1/3. At this parameter value a whole family of 
solutions exists describing 7r/2-front pairs with distances 
|a^2 — a^i I — 2|x| ranging from zero to infinity. In the non- 
gradient case these pair solutions propagate at speeds 



given by Eqn. (23). The degeneracy of solutions is lifted 
by higher order terms in Eqn. (0). 

Consider the gradient version of the amplitude equa- 
tion 



Bt=B + ]^B,,-\B\''B 



jiB*' + f^H{B,B*;d.,) 



(24) 



where H{B,B*;dx) includes higher order terms hke 
IBIS'S, \B\'^Bxx, etc.. The factor ^J. reflects the fact that 
fifth order terms in the amplitude equation are smaller 
by factor /i ^ 1 than the lower order terms. The effect 
of these terms is generally weak, but becomes important 
near 74 — 1/3. Consider, for example, the effect of the 
term S\B\^Bxx- Eqs. (0) include now the contributions 

^^^6{U^ + V^)Uxx and ^tiS{U^ + V^)Vxx , 

respectively. The corresponding contributions to 

Eqs. (H) and (|l|) are 

-^lSal [Bf^{x - xi) + Bl{x - X2)] B'f!,{x - xi) 



and 



-^5a2[Bl{x - xi) + BI{x - X2)\Bq{x - X2) 



The equation for the order parameter will now read 

27 Q 

X^~T7:dJ{x) + ^5^iK{x), (25) 



where 



/QO 
dz tanhz sech'^z tanh^(z + 2x) . (26) 
-00 



The integral (26) is elementary but the expression is 
lengthy and we do not display it here. The second term 
on the right hand side of Eqn. (Eq), whose origin is the 
fifth order term \B\'^Bxx, cannot be neglected in a /i- 
neighborhood of 74 ~ 1/3. Depending on the sign of 5 
two scenarios are possible as 74 is decreased. In both 
cases the x = (7r-front) solution is destabilized at 
74c = 1/3 + 8/^(5/21. When ^ > the x = solution 
is destabilized to a new pair of solutions, x±5 in a pitch- 
fork bifurcation. For |x| <C 1 the new solutions assume 
the approximate values, x± ~ ±-^^a/1 — d/dc, where 
dc = 74c — 1/3. When 74 is further decreased the two 
stable solutions x± move to ±00 on a 74 range of order 
ji. When (5 < bistability of the x — ^ solution and 
the X = ioo solutions first develop. As 74 is further 
decreased the x = solution becomes metastable until 
it loses completely its stability at 74c. Fig. || shows the 
potential 



V ^ 



9 



-dJ{z) - 6nK{z) 



dz , 



(27) 



associated with Eqn. (25) for both scenarios. 

The two scenarios are related by the symmetry d —f 
~d, S — > —5, t ^> —t oi Eqn. (pS]). The first scenario 
{S > 0) amounts to a pitchfork bifurcation from a stable 
X = solution to a pair of stable x± solutions that move 
to infinity as 74 is decreased. The second scenario {6 < 
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FIG. 5. The effects of the higher order term S\B\'^Bxx on 
the decomposition instability within the 4:1 resonance. The 
potential (kit) deforms from a single well to a single barrier 
as 74 is decreased past 74c. In the intermediate range two 
scenarios are possible: (a) For 5 > 0, the X = solution loses 
stability in a pitchfork bifurcation at 74c to a pair of solu- 
tions that move to ±00. Parameters: S = 1.0, /i = 0.01, 
74 = 0.339,0.337,0.336,0.335,0.333. (b) For S < 0, the 
X = solution remains stable while the x = ±00 solutions 
acquire stability and loses stability only below 74c. Parame- 
ters: \S\ = 1.0, fi = 0.01, 74 = 0.334, 0.332, 0.331, 0.330, 0.328. 
In both scenarios the deformations from a single well to a sin- 
gle barrier occur within a small range of 74 of order /i ^ 1. 
For comparison, an equivalent figure for the degenerate case 
{S = 0) is shown in (c). The only intermediate form between 
a single well and a single barrier is a flat potential, V = 0, 
ocurring at 7 = 1/3. 
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FIG. 6. Decomposition if a tt front into three 7r/3 fronts in 
the 6:1 resonance band. Parameters in Eqn. (bq): 76 ~ 0-9, 
jj,4 — —1.0, /i6 ~ —1.0. All other parameters are zero. 



0) amounts to a backward pitchfork bifurcation from an 
unstable x = solution to a pair of unstable x± solutions 
that move to infinity as 74 is increased. 

The higher order term \B\'^Bxx, and similarly other 
high order terms, lift the degeneracy of the lower order 
system (J^) at 74 — 1/3. For S > and in a small 74 range 
of order /i near 1/3 the instability becomes similar to the 
NIB bifurcation in the 2:1 resonance. But apart from 
the behavior in this small parameter range the overall 
behavior does not change: a 7r-front decomposes into a 
pair of 7r/2-fronts as 74 is decreased. 



IV. n-FRONT INSTABILITIES IN HIGHER 
RESONANCES 

We have found numerical evidence for the existence 
of similar tt- front instabilities within the 6:1 and 8:1 
bands. These findings suggest the following generaliza- 
tion: within the 2n:l band (n > 1) a 7r-front may lose 
stability by decomposing into n 7r/n-fronts. Consider the 
equation 



Bf 



1 



-Br 



{l + ii^o)B 



fi4\B\^B + fie\B\^B + ^is\BfB 
jiB*^ + -feB*^ + -fsB*\ 



(28) 



The normal form equation up to seventh order contains 
many more terms. Our purpose here, however, is just to 
demonstrate the 7r-front instability for some parameter 
values pertaining to the 6:1 and 8:1 bands. Fig. || shows 
the decomposition in the complex B plane of a 7r-front 
within the 6:1 band into three 7r/3-fronts. Fig. M shows 
the decomposition of a tt- front within the 8:1 band into 
four 7r/4-fronts. Fig. || shows a space-time plot of the 
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FIG. 7. Decomposition if a tt front into four n/4 fronts in 
the 8:1 resonance band. Parameters in Eqn. (pq): 78 = 0.75, 
/i4 = —0.5, fie ~ —0.5, /i8 ~ —1.0. All other parameters are 
zero. 



decomposition instability within the 6:1 band. The ini- 
tial unstable 7r-front decomposes into three 7r/3-fronts, 
traveling to the left or to the right depending on ini- 
tial conditions. Along with this process two intermediate 
phase states appear between the original white and black 
phases. 



V. IMPLICATIONS ON PATTERN FORMATION 

The TT-front instability in the 4:1 band has a pro- 
nounced effect on patterns. Despite the coexistence of 
four uniform phase states and the stability of 7r/2-fronts, 
asymptotic four-phase patterns appear only below the 
TT-front instability point 74 = 1/3. The reason is the at- 
tractive interactions between 7r/2-fronts when 74 > 1/3 
and the collapse into 7r-fronts. Thus, for 74 > 1/3 two- 
phase patterns prevail. These patterns form standing 
waves since tt- fronts are stationary. For 74 < 1/3 the in- 
teraction between 7r/2-fronts is repulsive and four-phase 
patterns prevail. These patterns travel since 7r/2-fronts 
propagate. 

Fig. |9p shows a stably rotating four-phase spiral wave 
for 74 < 1/3. Figs. y6, c, d show the collapse of this spi- 
ral wave into a stationary two-phase pattern as 74 is in- 
creased past 1/3. The collapse begins at the spiral core 
where the 7r/2-front interactions are the strongest. As 
pairs of 7r/2-fronts attract and collapse into 7r-fronts, the 
core splits into two vertices that propagate away from 
each other leaving behind a two-phase pattern. 
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FIG. 8. Decomposition if a tt front into three 7r/3 fronts 
in the 6:1 resonance band. The figures show space-time plots 
of numerical solutions of Eqn. (ESI) with parameters 75 — 0.9, 
/14 — —1.0, /i6 ~ —1.0, vo = 0.1. All other parameters are 
zero. 
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FIG. 9. Numerical solution of a two-dimensional version 
of Eqn. (0) showing the collapse of a rotating four-phase spi- 
ral-wave into a stationary two-phase pattern. The left column 
is |ylj and the right column arg(j4) in the x — y plane, (a) The 
initial four- phase spiral wave (computed with 74 < 1/3). (b) 
The spiral core, a 4-point vertex, splits into two 3-point ver- 
tices connected by a 7r-front. (c) A two-phase pattern develops 
as the 3-point vertices further separate, (d) The final sta- 
tionary two-phase pattern. Parameters: 74 = 0.6, fo = 0.1, 
a = (3 = 0,x= [0,64], y= [0,64]. 
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